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Abstract 

We study the conformation and scaling properties of a self-avoiding 
fluid membrane, subject to an osmotic pressure by means of Monte 
Carlo simulations. Using finite size scaling methods in combination 
with a histogram reweighting techniques we find that the surface un- 
dergoes an abrupt conformational transition at a critical pressure p* , 
from low pressure deflated configurations with a branched polymer 
characteristics to a high pressure inflated phase, in agreement with 
previous findings g, |. The transition pressure p* scales with the 
system size as p* oc N~°^, with a = 0.69 it 0.01. Below p* the en- 
closed volume scales as ^ oc A^, in accordance with the self-avoiding 
branched polymer structure, and for p \ p* our data are consistent 
with the finite size scaling form V oc N^+, where /3+ = 1.43 it 0.04. 
Also the finite size scaling behavior of the radii of gyration and the 
compressibility moduli are obtained. Some of the observed exponents 
and the mechanism behind the conformational collapse are interpreted 
in terms of a Flory theory. 
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1 Introduction 

Over the past several years there has been an increased interest in the phase 
behavior and morphological properties of flexible, fluid interfaces. Beside 
the theoretical challenge in understanding surfaces as two dimensional gen- 
eralizations of polymers, they are expected to be of relevance in physical 
systems ranging from simple interfaces between coexisting fluids close to 
the consolute point to surfactant interfaces in e.g. microemulsions and lipid 
membranes 0, |^, H, 0]. 

In the context of modelling lipid vesicles it is of particular importance to 
study the conformation of closed membranes in the presence of surface ten- 
sion, bending rigidity, osmotic pressure between the interior and the exterior 
of the membrane, etc. For a simple, closed, fluid membrane with bending 
rigidity k ■ ksT, pressure p ■ ksT and fixed overall topology, the vesicle con- 
formation is determined by the energy functional, 



where Ri and R2 are the principal radii of curvature P]. This phenomeno- 
logical form is the basis for detailed descriptions of shape transformations 



ciently large, unconstrained fluid self-avoiding surfaces collapse into branched 
polymer-like structures with the characteristics of branched polymers for 
T > 0|jTl[. This general feature will not be changed by e.g., the presence 
of bending stiffness. However, the osmotic pressure difference between the 
interior and exterior of a closed surface is a quantity of direct experimental 




(1) 



of rigid vesicles with k ^ 1 |^ 10 1. There is a general consensus that suffi- 
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relevance, which is expected to change this situation. 

The properties of two-dimensional vesicles (ring-polymers), subject to a 
pressure p controlling the enclosed volume (area), have been explored in a 
number of studies [|12|-[^. Leibler, Singh and Fisher |T2| verified by nu- 
merical studies that 2-D flexible vesicles change continuously from a deflated 
state with characteristics of branched polymers to an inflated state as p is 
increased. This change is accompanied by crossovers in scaling exponents 
characterizing the surface. The scaling-behavior of the inflated regime can 
be interpreted in terms of Pincus stretching exponents for polymers pO|] . 
Real-space-RG studies of a lattice model equivalent to the 2-D vesicle prob- 
lem |^0| and exact configurational enumeration techniques [T^ show agree- 



ment with the numerical studies. The universal relations have also been 



confirmed by use of conformal techniques on the lattice model |18|. Further- 
more, analytical studies of pressurized, flexible, but self-intersecting surfaces 
in 2-D have been carried out Only few studies have been carried out 



in 3-D. In |^ ^ pressurized, flexible, fluid surfaces have been considered 



in the context of the confined phases of some lattice gauge models. These 
studies have focussed on the effect of pressure on vesicle topology for deflated 
vesicles. 

It has recently been demonstrated by numerical studies that pressurized, 
flexible surfaces in 3-D undergo a first-order-like transition from a deflated 
phase with branched polymer behavior to inflated configurations. Gompper 
and KroU based their studies on a triangulated random surface model [Q], 
while Baumgartner apphed a plaquette model for a self-avoiding surface P, §], 
where the surface is composed of the domain walls of a spin-model. In these 



studies the basic scaling relations, established for 2-D flexible, pressurized 
surfaces, were generalized to three dimensions. 

In the present paper we also consider a fluid membrane with an osmotic 
pressure p, but with zero bending rigidity n. Its properties are analyzed on 
the basis of computer simulations of a triangulated random surface represen- 
tation of the fluid membrane. Using the Ferrenberg-Swendsen reweighting 
techniques we attempt a precise determination of the critical pressure 
and the scaling exponents in the deflated and inflated phases. 

2 Model 

The model considered here is a simple extension of the Ho-Baumgartner de- 
scription of self-avoiding surfaces It consists of a closed triangular net- 
work of beads positioned at the vertices of the network with the topology 
of a sphere. To each vertex position Xj is associated a bead with hard core 
diameter 1. In order to enforce the self- avoidance constraint the length of the 
flexible tether between neighboring vertices is chosen to be less than ^/2. The 
fluidity of the surface is modelled by a dynamical change of the connectiv- 
ity between the vertices by means of 'flip' transformations which correspond 
to deleting a tether or link between neighboring vertices and attempting to 
form a new link between the adjacent vertices in the two triangles involved. 
Finally, the shape changes are generated by a local updating of the vertex 
position, the 'shift' transformation, Xj Xj + 5Xj, where 5Xj represent 
an incremental change in the local surface position. 5Xj is picked randomly 
within a cube, where the cube-side is adjusted to preserve self-avoidance. 
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These procedures ensure sampling over all possible self-avoiding, piece-wise 
linear surfaces consisting of beads and with the topology of a sphere. The 
updating of the surface shapes is posed by standard Monte Carlo techniques 
on the partition function 

Z= ^M-H/kBT). (2) 

Surface shapes 

In the simulations a discretized version of the Hamiltonian from Eq.(|l|) 
with K = is applied. In order to obtain a fast updating rate a dynamic 
sublattice structure is implemented which keeps track of the N vertices during 
shape transformations. We furthermore use data structures consisting of 
linked lists which, given an arbitrary vertex or link, allow an identification of 
the adjacent triangles. The last feature implies an oriented triangulation and 
leads to the following expressions for the area of the surface and the enclosed 
volume: 

A = Y.^A, = W I {X, - X,) X {Xi - X,) I, (3) 

n ^ ijl 

and 

V = J2AV^ = Wx^■iX,xXl). (4) 

The indices ijl pertain to the corners of an oriented triangle and the sum 
runs over all triangles of the triangulation. 

Using a stress ensemble we have collected data for a range of p values and 
values (A^ is the number of vertices). In a simulation we typically measure 
the following quantities: 

• We monitor the shift and flip rates in order to ensure a suitable balance 
between fluidity (the flip rate) and shape changes (the shift rate) 
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• For the characterization of the over-all membrane size we sample the 
radius of gyration Rq: 

1 ^ - 

Rl^^Y.i^^-XcM)^ (5) 
i=l 

where Xcm — jf J2iLi is the center of mass. 

• the area A and the volume V of the membrane are monitored. 

The measured quantities are stored for construction of probability distri- 
butions and evaluation of equilibrium thermal averages. It is characteristic of 
fluctuating membrane systems that the correlations times are quite long and 
grow fast with the system size. Careful analysis of the time series is thus very 
important to ensure proper equilibrium sampling of the fluctuating surface. 
E.g., for = 400 we sample 20 Mill MCS/vertex in order to equilibrate the 
system and perform measurements over about 400 Mill MCS/vertex. We 
have investigated system sizes ranging from N — 4:9 to N — 400. 



3 Simulations and Data Analysis 

The Ferrenberg-Swcndsen reweighting of the probability distributions allow 
for an approximate determination of thermal averages for a pressure value p 
close to the simulation pressure po, if the phase space is properly sampled, 
e.g. for the averaged volume (V) it takes the form: 

_ j:iViP{Vi)exp{-ApVi) 

where Ap = p—po, and PiVi) is the sampled distribution. The sum is running 
over all measured volume values. Figure 1 shows the averaged vesicle volume 
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versus pressure p for two different system system sizes = 100 and 196, ob- 
tained from simulations at different pressures. The lines are the results of the 
reweighting technique, using only one simulation near the assumed critical 
pressure. The dramatic change in the vesicle volume over a narrow pressure 
interval indicates the presence of a phase transition in the thermodynamic 
limit N ^ oo from a low pressure deflated phase to a high pressure inflated 
phase in accordance with the simulations of Gompper and Kroll The plots 
of volume versus pressure thus allow an approximate determination of the 
critical pressure p*. The reweighting technique also enables us to reconstruct 
the phase diagram. In the insert of Fig. 1 we show the reconstructed volume 
versus pressure curves for = 81,100,121,144,169,196,225,256,324,400 
obtained by use of Eq.@. 

A more precise determination of the critical pressure can be obtained by 
evaluating the total volume fluctuations {V"^) — {V)"^ and estimate the peak 
position. Figure 2 shows the measured p* for different system sizes, p* is 
well described by a power-law in A^: 

p* ~ (7) 

where a = 0.69 ± 0.01, obtained from a linear fit to the double-logarithmic 
representation of the data. All indicated error bars are related to the regres- 
sion analysis. 

The Ferrenberg-Swendsen reweighting technique can also be applied 
for the construction of the probability distribution P{V) at specific pressures 
near the anticipated critical pressure p*. In Fig. 3 we show examples of mea- 
sured histograms representing P{V) and P{Rq) for A^ = 144 at p = 0.52 and 
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p = 0.54. The measured P{V) at p = 0.52 is in good agreement with the 
probabihty distribution obtained from reweighting the histogram measured 
at p = 0.54. At p* P(y) has a characteristic double-peak structure repre- 
senting the two phases. Lee and Kosterhtz used this technique to determine 
finite-size transition points for lattice models by matching the peak-heights, 
and by investigation of their finite-size behavior information about the na- 



ture of the transition in the thermodynamic limit can be obtained I^Bj. It is 
a key assumption in this procedure that the probability distribution for some 
density can be well described as the superposition of Gaussian distributions. 
This is not fulfilled for PiV), but some properties of the distribution will be 
considered. P{V) can in a first approximation be written as the superposition 
of contributions from infiated and defiated contributions: 

P{V) ~ a_P_(\/) + a+P+{V), (8) 

as suggested by e.g., Fig. 3. The maximum probability volumes for the two 
peaks V- and V+ can easily be estimated and is shown in Fig. 4. We find 
the dependence: 

V±~iV^±, (9) 

where the associated scaling exponents are /3_ = 0.99 ± 0.01 and /?+ = 
1.43 ±0.04. 

Since the two peaks in P{V) display very different dependences of A^, P_ 
and P+ have been analyzed separately. In Fig. 5 P+ is given in a double log- 
arithmic representation up to a multiplicative factor. The plot indicates that 
the volume distribution for the inflated phase of the surface has a universal 
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form for large N. 

PAV) j^M^h (10) 

where /+ is represented in Fig. 5 as the hmiting distribution as becomes 
large. Our data are not sufficient for a thorough analysis of the form of /+. 

The probability distribution describing the volume fluctuations of the 
deflated configurations are given in the insert of Fig. 5. Similar to P+, P_ 
approach a universal form as N becomes large: 

p_(y)~/_(y-y_), (ii) 

/_ is represented in the insert. The volume fiuctuations ((Ay)^)_ are thus 
independent of N. 

At p* the probability distribution for the radius of gyration Rq displays a 
double-peak similar to the probability distribution of the volume, as shown 
in Fig. 3. The peak-positions are identified with the radii of gyration of the 
infiated and the defiated phases R"^. The calculated R"^ shown in Fig. 6 
indicate a relationship: 

i?|^Ar2^±, (12) 

where 2i/_ = 1.05 ± 0.05 and 2i/+ = 0.81 ± 0.01. 

In Fig. 7 the data for the total volume fiuctuations K — {V^) — {V)'^ at 
p* is plotted against N. From the regression analysis in the log-log represen- 
tation we obtain a dependence: 

K{p*)^N\, (13) 

where 7 = 3.62 ± 0.02. 
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4 Discussion and Conclusion 



The above analysis confirms that fluid, flexible vesicles undergo an abrupt 
conformational change between a deflated and an inflated regime at some 
small pressure value which diminishes with the system size according to a 
power-law, p* oc A^~", where a = 0.69 ± 0.01. The value of this exponent is 
in disagreement with previous flndings 0.5 and 1 0. 

In the deflated regime p < p* the data are consistent with branched 
polymer conflgurations of the surface: 

(X {R'g)^ (X AT, (14) 

and thus support the general expectation that unpressurized, fluid membrane 
conformations are controlled by the branched polymer flxed point, i.e., the 
stable K = crumpling flxed point p9| . 

In the inflated phase for p > p*, generalizing Pincus' result for polymer 
stretching Gompper and Kroll fl]] and Baumgartner P| propose the 
scaling form 

{V)+ ~ p3{2P-l)^3P^ (^5) 

where z/ is a new scaling exponent characterizing the stretched phase. In the 
following we will use the short notation R = ^ {Rq) + . 

At the critical pressure inserting the scaling result p* ~ A^~" in Eq. ([T5|) , 
we obtain for {V)+ in the inflated phase, 

~ Ar3[^-"(2--i)l, (16) 
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and we derive the relationship (3+ = 3[z/ — a{2u — 1)]. From the numerically 
obtained exponents /?+ ~ 1.43 and a ^ 0.69 in Eqs.(|) and (^) we find 
u ~ 0.56 which is in a reasonable agreement with Gompper and Kroll 
{u = 7/12 = 0.583). However, a similar analysis of Rq leads to a significantly 
different estimate u ~ 0.75. 

Another consequence of Eq.(|l5D is, that for p \p*: 

which is contradictory to the numerical finding ((AV)^)+/(1^)+ ~ 
where ((Al^)^)_|_ ~ A^2/3+ ^j^g transition is a consequence of Eq. (p!0|) . Fi- 
nally, we find {V)j^/ B? oc A^^+~3'^+^ where /3+ — 3z/+ ^ 0.22 at the transition, 
so the relation (y)+ R? is not consistent with our data. The numerical 
data, obtained in the transition region, are not sufficient for an extended 
scaling analysis, however, with the above considerations we conclude that 
they do not support an ansatz of the type Eq. ([T5|) . 



In the following we will discuss the infiation-defiation collapse transition 
of a fiuid surface within the framework of a Flory theory which has been 
applied to the description of polymeric, fiexible surfaces and attempted 



applied in the context of plaquette surfaces [Q. The application of sufficient 
pressure will restrict the configurational phase space and prevent branching 
of the surface. Such a surface can for "small" pressures allow an effective 
description in terms of a Flory theory for the free energy: 

F. = - + t;-^-p7iVi?. (18) 
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The first two terms represent the standard Flory theory for a surface, where 
Ro is the radius of gyration for a self-intersecting surface, Ro ~ N'^°, where 
z/q = 0, and v is the excluded volume parameter, d = 3 is the dimension of 
the embedding space. In the last term the pressure couples to the volume: 

{V)+(x{f dAn-{X- Xcm))+ oc (cos 9)+ ■ NR, (19) 

J A 

where 9 = 1{X — Xqm, and n is the the surface normal at surface position 
X. Here it is used that \X — Xcm| is distributed around R in the inflated 
configurations, so Eq.(|l9]) represents the leading contribution to {V)+. The 
parameter 7 in Eq.(0) is thus proportional to (cos 6*)+ > 0. The stationary 
condition for Fj at p = can be characterized by a critical exponent u: 

R (X N"", u= = - = 0.4, 

d + 2 5 

(V^)+ oc NRocNf^, (3 =l + v = - = 1.4. (20) 

5 

These Flory exponents are indeed close to the numerically obtained v+ = 0.41 
and /5+ = 1.43. But since Fi oc N'^'^, these configurations are unstable with 
respect to branched polymer collapse for large and must be considered 
as a metastable free-energy branch for p = 0. The branched polymer config- 
urations are completely entropy dominated with an approximate free energy 
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Ffep = - \niz*)N + ^ In(iV) - pN, (21) 

where z* > 1. Let us now apply a small pressure, so that the relations in 
Eqs.(pnD still hold. The characteristic pressure p* separating the unstable 
p = region from the inflated region in Fi is given by p*N^^'^ ~ N'^'^ leading 
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to p* oc A^~", where a = 1 — u = 3/5 = 0.6. For a choice of a smaller 
than 3/5 the perturbative considerations break down and the pressure term 
will dominate in Eq.(|l^), which eventually leads to complete inflation with 
F oc —pN^f"^ for ^ oo. For a larger than 3/5 the surface will be controlled 
by the p = behavior as becomes large. This simple analysis thus sug- 
gests that the trigger of the inflation-deflation transition can be understood 
on basis of the cross-over from the flaccid to the inflated conformations of 
closed, flexible, polymeric surfaces under pressure. The strong fluctuations in 
the transition region, e.g., observed in the probability distributions, are not 
included in this consideration. In particular it must be expected that fluctua- 
tions in the weakly inflated regime will modify the conditions for the stability 
of the inflated configurations. The deviation of a from the numerically ob- 
tained exponent a = 0.69 is thus not surprising in light of the crudeness of 
the approximations involved. 

In the present paper we have reanalyzed the deflation-inflation transition 
of a fluid membrane subject to an inflating pressure discovered by Gompper 
and Kroll [||], and Baumgartner p|. In the transition region we find low- 
pressure configurations with the characteristics of branched polymers and an 
exponent for the radius of gyration z/_ = ubp = 0.52±0.03 and high pressure 
configurations, with a characteristic exponent z/+ = 0.41 ± 0.01. Analysis of 
the volume distribution functions show that the two phases display universal 
properties at p*, characterized by exponents for the volume /?_ = 0.99 ±0.01 
and /?+ = 1.43 ± 0.04 for the two types of configurations in the transition 
region. Further, a simple analysis of the transition within the framework 
of Flory theory accounts for some of the observed exponents and suggests 
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a mechanism for the transition. The exponents characterizing the finite-size 
behavior of the transition pressure and the total volume fluctuations are sub- 
jects to further investigations. 
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Figures 



1. In this figure the volume averages (V) obtained from simulations at 
different pressures are plotted vs p, for two different system sizes (A^ = 
100 and N — 196). The lines are the results from the Ferrenberg- 
Swendsen reweighting for a single simulation at p = 0.7 for N — 100 and 
p = 0.445 for N = 196. The insert shows the volume values obtained 
with the Ferrenberg-Swendsen reweighting technique from simulations 
near the critical pressure, for N = 81, 100, 121, 144, 169, 196, 225, 256, 
324, 400. 

2. This figure shows log(p*) vs log(A'"). The values forp* are obtained from 
the peak in the total compressibihty ((AF)^), which was calculated 
using the Ferrenberg-Swendsen reweighting technique. The exponent 
is obtained from fitting to —0.69 ± 0.01. 

3. The distributions P{V) and P{Rg) are plotted for two different pres- 
sure values and the system size = 144. The sharp peak in the P{Rg) 
distribution belongs to the inflated phase, whereas the corresponding 
peak in the P{V) distribution is broader. For the branched polymer 
configurations it is opposite. 

4. The volumes V_ and V+ corresponding to the two peaks in P{V) at 
p* are plotted against in a double logaritmic representation. The 
regression fines correspond to V± oc N^"^, where /3_ = 0.99 ± 0.01 and 
(3+ = 1.43 ± 0.04. 
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5. The unnormalized volume probability distrubutions at p* for N = 100, 
121, 144, 169, 196, 225, 256, 324, 400. P{V) is for each N multiplied 
by a common factor so the maxima are equal unity. P+: The part 
of the histograms associated with the inflated configurations is plotted 
against V/N^+ in a log-log representation. P_: In the insert is the low- 
volume parts of the histograms representing the inflated configurations 
similarly given in a semi- logarithmic plot. 

6. The radii of gyration and R\ at p* are plottet against A^. The data 
are indicative of a relationship oc A^^± where — 1.05 ± 0.05 and 
i/+ = 0.81 ± 0.01 

7. The maximum of the total volume fiuctuations {{/^Vf) = {V^) - {Vf 
(corresponding to p*) is plotted vs. in a log-log representation. The 
fit indicates ((AT/)^)^^^ ~ N'^ , with 7 = 3.62 ± 0.02 
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